Chemical or biological activity in open chaotic flows 
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We investigate the evolution of particle ensembles in open chaotic hydrodynamical flows. 
Active processes of the type A + B — > 2B and A + B —* 2C are considered in the limit of weak 
diffusion. As an illustrative advection dynamics we consider a model of the von Karman vortex 
street, a time periodic two-dimensional flow of a viscous fluid around a cylinder. We show that 
a fractal unstable manifold acts as a catalyst for the process, and the products cover fattened-up 
copies of this manifold. This may account for the observed filamental intensification of activity 
in environmental flows. The reaction equations valid in the wake are derived either in the 
form of dissipative maps or differential equations depending on the regime under consideration. 
They contain terms that are not present in the traditional reaction equations of the same active 
process: the decay of the products is slower while the productivity is much faster than in 
homogeneous flows. Both effects appear as a consequence of underlying fractal structures. In 
the long time limit, the system locks itself in a dynamic equilibrium state synchronized to the 
flow for both types of reactions. For particles of finite size an emptying transition might also 
occur leading to no products left in the wake. 

(February 5, 2008) 



I. INTRODUCTION 

Active processes taking place in chaotic hydrodynam- 
ical flows have attracted recent interest fj]]- (h]]. By 
chaotic we mean time-dependent but nonturbulent ve- 
locity fields with chaotic tracer dynamics (Lagrangian 
chaos) O]- |L3|]. In the simplest approximation we can 
assume that the advected particles undergo certain chem- 
ical or biological changes but do not modify the fluid flow. 
The motivation for such studies has been to understand 
the effects of imperfect mixing [B due to the underly- 
ing chaotic particle dynamics. The implications can be 
perceived in laboratory experiments p|,|l4|], but the ef- 
fects are perhaps more striking in environmental flows. 
In particular, there is increasing evidence of filamental 
structures in the product distribution of environmental 
processes both in the atmosphere, such as ozone reac- 
tions |§-[||]U , and resulting from the interaction of pop- 
ulations of microorganisms in the sea, such as plankton 
distributions Jl5|,[l6| ■ Our aim is to show that these struc- 
tures might be consequences of the fractal structures of 
the reaction free flows j7j. 

Here we shall consider open flows with asymptotic sim- 
plicity in which the velocity field in the far up and down- 
stream regions is uniform. A well-known (time-periodic) 
laboratory example is the flow around a cylinder. Its ac- 
tual realization can be observed in environmental flows, 
like, e.g., in the fluid motion in the wake of a pillar or 



in the motion of air behind an isolated mountain. A 
unique feature of such open flows is the pronounced and 
stable fractal feature associated with the chaotic tracer 
dynami cs JlTf ]- p7| , which is clearly measurable in exper- 
iments pq] . The central object governing the dynamics 
is a nonattracting chaotic saddle [ p9| containing an in- 
finite number of periodic and nonperiodic tracer orbits 
which remain bounded and never reach either the far 
up or the downstream regions. A characteristic quantity 
of the saddle is its escape rate k whose reciprocal value 
is the average chaotic lifetime. The far up and down- 
stream regions are foliated by the saddle's stable and 
unstable manifold, respectively. The saddle's unstable 
manifold directs tracers ever approaching the saddle to 
the far downstream region. Though both the saddle and 
its manifolds are not space-filling fractal objects, it has 
been pointed out |l7j - j28j that the unstable manifold is 
the avenue of propagation and transport in such flows. It 
is the pronounced fractal structure of such flows (which 
is not present in closed flows) that makes them specially 
interesting catalyst of active processes |7| . 

In this paper we consider the advection of active par- 
ticles in flows with asymptotic simplicity in which the 
activity is assumed to be of chemical or biological origin 
in the simplest possible form. The reaction is a kind of 
"infection" leading to a change of certain properties, such 
as color of reacting particles. Particles with new prop- 
erties are the products. Since it is in the close vicinity 
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of the chaotic saddle and its unstable manifold that the 
particles spend the longest time close to each other, it is 
there where the effect of the activity is most pronounced. 
It is then natural to expect that the products should ac- 
cumulate along the unstable manifold and trace out this 
fractal object. 

In our work we support this conjecture and present 
a detailed analysis of such active processes. We show 
that the unstable manifold of the chaotic saddle is the 
backbone of the reaction. The newly born components 
cover the branches of the unstable manifold with a well 
defined average width e* . Thus, an effective fattening up 
of the fractal set takes place due to the activity of the 
tracers. This implies that on linear scales smaller than 
this width e*, fractality is washed out, but a clear frac- 
tal scaling of the material with a dimension Dq can be 
observed on larger scales. This fractal dimension is the 
same Do as that of the unstable manifold in the reaction 
free flow. Although the fractal set itself is a set of mea- 
sure zero, the amount of chemical products is finite due 
to the fattening-up process of this manifold. 

A consequence of the fractal backbone is that the 
amount of the reaction product follows a singular scaling 
law with irrational Do-dependent powers of the number 
of product particles, signaling a singular enhancement of 
productivity . (The enhancement of activity is meant 
in comparison with non-chaotic, e.g., stationary flows.) 
This singularly enhanced rate of activity has profound 
practical consequences. It may account for the observed 
filamental patterns of intense activity in environmental 
flows an effect that cannot be explained if one 

considers diffusion processes alone. In this work we show 
how small-scale structures are generated in the dynamics 
of active particles, and how these dynamical structures 
are responsible for the enhancement of activity. 

In summary, the effect of the chaotic saddle producing 
this activity is twofold: (i) to keep the reacting particles 
longer, as given by the escape rate k, in the interaction 
region, and (ii) to generate the unstable manifold's fractal 
structure that brings the reacting particles together. 

We derive the corresponding reaction equations in the 
form of maps or differential equations depending on the 
regime under consideration. Such processes are general- 
ization of classical surface reactions pOj, but, by contrast, 
in our case the surface is a fattened-up fractal. The re- 
action equations contain new terms not present in the 
traditional well-stirred reaction model of the same pro- 
cess. In spite of the passive tracers' Hamiltonian dynam- 
ics, these reaction equations turn out to be of dissipative 
character possessing attractors. 

We find that the chemical activity and the advection 
by the hydrodynamical flow are in permanent competi- 
tion. Due to this competition, most typically, a kind of 
stationary state sets in after sufficiently long times. In the 
case of time-periodic flows of period T, the asymptotic 
state is typically also periodic with T, i.e., the reaction 
becomes synchronized to the flow. Thus, a chaotic par- 
ticle dynamics is consistent with a non-chaotic reaction 



dynamics. 

To be more specific, we consider simple kinetic mod- 
els Q] with disk-like particles. Two particles of differ- 
ent kinds undergo a reaction if and only if they come 
within a distance cr, which is the interaction range. Due 
to the incompressibility of the fluid (which is always a 
good approximation for velocities much below the speed 
of sound), two-dimensional flows are area-preserving. 
Therefore, if we cover a layer of fluid by particles, this 
single-particle coverage of the layer will not change in 
time. In other words, the average linear distance Eq < a 
between the actual nearest- neighbor particles is constant. 
(Note, however, that the linear distance between specific 
particles is typically increases in time due to the chaotic 
dynamics.) Thus, the result of the reaction is just the 
spreading of the products in the layer of particles being 
advected by the flow. We emphasize again that parti- 
cles are assumed to have no feedback on the flow. Fur- 
thermore, the advection dynamics is purely determinis- 
tic, i.e., we work in the limit of weak diffusion where the 
reaction range represents some kind of diffusion distance 
too. 

We shall consider both an autocatalytic process, A + 
B — > 2B, and a collisional reaction A + B — > 2C. In 
both cases A is considered to be the background material 
which covers initially the full infinite layer of observation. 
In the autocatalytic process a single seed of particle B 
is sufficient to trigger reactions leading to the station- 
ary state, while in the collisional reaction a continuous 
feeding of material B is necessary. 

For computational convenience we assume that the re- 
actions are instantaneous and take place at integer mul- 
tiples of a time lag r. We shall see that an important di- 
mensionless parameter will be the ratio between the time 
lag and the average chaotic lifetime: v = r/(l//t) = tk. 
This can also be considered to be the dimcnsionless re- 
action time, whose reciprocal value tells us how many 
reaction events occur on the characteristic time of chaos. 

The case of time continuous reactions is obtained in 
the limit r — > (or v — > 0) by keeping the reaction front 
velocity v r finite. In this limit we assume, too, that the 
average distance Eq between particles goes also to zero, 
and we obtain a continuous distribution of particles. We 
call this limit the chemical frame. A fractal product dis- 
tribution is then expected to appear if the reaction is slow 
compared to the flow (v r <C kL with L as a characteris- 
tic length). An example for time continuous reactions is 
related to the depletion of ozone at the polar vortex: the 
trimolccular reaction of CIO with NO2. In late winter 
and early spring the polar vortex exhibits high concen- 
trations of CIO and very low concentrations of NO2 while 
outside the vortex the situation is typically reversed with 
relatively high concentrations of NO2 and low concentra- 
tions of CIO. Thus the reaction C10+N0 2 -> C10N0 2 is 
a natural candidate to produce a filamental CIONO2 dis- 
tribution along the edge of the polar vortex ||,||,|9|,[5l],[52| 
on the time scale of a few days, where the molecular dif- 
fusivity is negligible. 



2 



In order to contrast our work with the conventional 
description of chemical processes, we briefly discuss now 
the continuous time dynamics of the autocatalytic sur- 
face reaction A + B — > 2B in a uniform flow. Let us 
observe the flow moving to the right with velocity vq in 
a unit square (fixed to the observer at rest). A seed par- 
ticle of type B is kept fixed about the middle of the left 
boundary. Particles of type A are distributed with uni- 
form density everywhere on the surface of the flow (also 
upstream). The seed particle starts to interact with its 
A neighbors transforming them into B. Since B par- 
ticles are transported away, more and more A particles 
are converted into B. Let us assume that at time t the 
B particles cover a triangle across the square which is 
symmetric about a horizontal line (Fig. ^). 




FIG. 1. Schematic diagram of a uniform flow of velocity vo 
with autocatalytic reaction. A single seed B is kept fixed at 
the left corner of the gray triangle occupied by material B 
which lies in a layer of white A particles. The observation 
region is a unit square (dashed line). 

We assume that the half angle a is small. The area 
Ab occupied by B is simply a. The change in this area 
during time dt is due to a horizontal displacement v dt 
and a vertical increase v r dt of both fronts, where v r is 
the reaction front velocity. The gain of the area Ab is 
just 2(— voa + v r )dt = 2(— vqAb + v r )dt. The differential 
equation governing this area is thus 

A B = -2v Q A B + 2v r . (1) 

It has a stationary solution, corresponding to a stabilized 
triangular distribution of B particles of area A* B = a* = 
v r /vQ. We shall see that the presence of a saddle in a time 
dependent nonuniform flow results in a slower decay and 
a faster production (2vq is replaced by the escape rate 
n and the production term will contain a factor with a 
negative power of the area itself due to the fractality of 
the unstable manifold). 

For generality, we also investigate cases where the time 
lag t is finite so that its dimensionless version v is of order 
unity or larger. If v exceeds a critical value, we find that 
no product remains, i.e., an emptying transition takes 
place. Because of the finite values of r, and the discrete 



character of the particles (eq might also be considered as 
the size of particles), this latter effect might be of rel- 
evance to biological processes accompanying advection. 
An example can be a crude model of the dynamics of 
plankton populations |^,[^] in the presence of a time 
dependent flow. The so-called zooplanktons (B) have a 
daily rhythm: they sink down during night time but come 
up to the surface of the sea again during day-time when 
they eat up phytoplanktons (A) 7 reproduce themselves 
and then grow in number. 

This paper is organized as follows. In Sec. II, we 
present the model and the numerical procedure used, 
while the results for both reaction types are shown in 
Sec. III. A detailed theory based on these observations is 
derived in Sec. IV. The concluding Sec. V gives remarks 
on properties expected to be valid in more general models 
of active processes in open flows. 

II. THE MODEL FLOW AND NUMERICAL 
PROCEDURE 

The flow chosen to illustrate the fractal active dynam- 
ics is an example of a two-dimensional, incompressible 
time-periodic fluid motion, the case of the von Karman 
vortex street in the wake of a cylinder [psft . The 

radius R of the cylinder and the period T of the flow are 
taken as the length unit and the time unit, respectively. 
In what follows we keep the flow parameters constant, im- 
plying a fixed value of the escape rate k, and investigate 
the dependence of the reaction outcome on parameters 
like the reaction range (a) and time (v). 

For simplicity, we use an analytic model for the stream 
function introduced in Ref. |2(J (the explicit form of the 
stream function can also be found in the Sec. Ill of Ref. 
p4j). This model has been motivated by direct numer- 
ical simulations at Reynolds number of about 250 [fl9|| , 
and has been used successfully to reproduce qualitative 
features of the tracer dynamics. The escape rate of the 
particles in the reaction free flow is k — 0.36 and the 
fractal dimension of the unstable manifold is Dq = 1.61, 
while the background flow velocity is vq = 14 p(i| , p4| . 

Since the flow is periodic, we fix the "phase" of the re- 
action relative to the flow. We consider time zero, t = 0, 
to be the instant when a vortex is born close to the sur- 
face of the first quadrant of the cylinder and, simultane- 
ously, a fully developed vortex is detaching in the fourth 
quadrant of the surface |20[|24|| . 

For convenience, we carry out the simulations on a 
uniform rectangular grid of lattice size Eq covering both 
the incoming flow and the mixing region in the wake of 
the cylinder. This £o also corresponds to the average 
distance between nearest-neighbor particles. If there is 
a tracer inside a cell, it is always considered to be in 
its center; the reaction range is bounded from below by 
the lattice size: £o < c. This projection of the tracer 
dynamics on a grid essentially defines a mapping among 
the cells. 
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The course of the chemical reaction starts with nearly 
all cells occupied by species A, the background material. 
Few cells contain B distributed according to the initial 
conditions chosen for the type of reaction under consider- 
ation. One iteration of the process just described consists 
of two mappings in involution. The first mapping models 
the advection of the particles on the chosen grid, while 
the second models the instantaneous active process (e.g. 
chemical reaction) occurring on the same grid of cells. 
In fact, in any closed region considered there is a loss of 
the products due to the advection but also a gain in the 
product amount due to the reaction. Later we find a bal- 
ance between these two competing effects manifested in a 
steady state in the product distribution. The simulation 
consists of a repeated application of advection and re- 
action steps. We apply different algorithms for different 
reaction types. 



A. Autocatalytic reaction: A + B — > 2B 

If a tracer starting from the center of a cell is advected 
into another one after time t, then the latter cell is con- 
sidered to be the image of the first one with respect to 
the dynamics. After an application of the map, a cell will 
be considered occupied by reagent B if it is an image of 
at least one B cell. Otherwise the cell is considered to 
contain species A after the mapping. In addition, if a cell 
contains B at the time of the reaction, all of the 8 neigh- 
boring cells are infected by B. Consequently, Eq plays 
the role of the interaction range a in our simulation. 



B. Collisional reaction: A + B — > 2C 

In this cell is considered to be the image of 

another one with respect to the dynamics if its center's 
preimage is inside this other cell r time earlier . This 
defines the mapping among the cells due to the dynamics. 
After the action of the mapping, the reaction can modify 
the cell contents: Any cell containing A (B) before the re- 
action becomes C, if there is a B (A) cell within a radius 
a from its center. Otherwise the cell keeps its content. 
Numerically we found convenient to store the configura- 
tion of the lattice just before reactions only. Then the 
content of a cell at the time just before a reaction can 
be deduced from its preimage and the neighbors of the 
preimage according to the following criterion: If, r time 
earlier, among the preimage cell and its neighbors there 
were both types A and B present, then the cell must 
have become C during the last reaction; if all of them 
were of one type only (apart from C, which is inert), the 
cell inherits the type of its preimage. This means that 
we unify the advection-reaction process in one mapping 
connecting the cell contents just before reactions. In all 
experiments the reaction range a is on the order of the 
lattice size Eq. 



III. RESULTS 



A. Autocatalytic reaction: A + B 



2B 



Initially, we introduce a seed of reagent B in front of 
the cylinder. Since there are only two species in the 
system, we monitor only reagent B. Values referring to 
material A inside the computational domain can be ob- 
tained from mass (in our two-dimensional model, area) 
conservation. Figure |^ displays the spreading of reagent 
B (black) in the course of time. Note the rapid increase 
of the B area and the quick formation of a filamcntal 
structure that becomes stationary after a few time units, 
but changes periodically with the period of the flow. 
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FIG. 2. Time evolution of a seed of 400 B particles (black) 
placed in the flow in front of the cylinder (on a square of 
linear size 0.2) at time 0. As the auto-catalytic reaction 
evolves having A as the background particles (white), the 
amount of B increases and traces out a complicated object 
in the wake of the cylinder. After some initial increase, a 
steady state sets in. The snapshots (a-j) are taken at times 
t = 0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.6, 2.0, and 3.0 respectively, 
right before a reaction takes place. The computational do- 
main —3 < x < 5 and —1.5 < y < 1.5 covers both the in- 
coming and mixing regions. The time lag between reactions 
is r = 0.2, consequently v — 0.072, and the lattice size is 

£ = cr = 0.01 

To support this qualitative observation, Fig. [| shows 
the number of B particles in the computational domain 
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as a function of time. After four periods, a self-repeating 
time dependence sets in. 




FIG. 3. The dependence of the number Ag (t)/eq 
(eo = a) of B particles in the computational domain of Fig. 2 
on time t — nr right before the reaction events. Note the sta- 
tionary time-periodic behavior reached after about four time 
units. 

This means that the chemical reaction takes over the 
flow's basic periodicity and reaches a stationary state: 
the number of cells being born in the reaction is the same 
as the number of cells escaping due to the advection dy- 
namics. 
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FIG. 4. (a) Time evolution of a droplet of the same size as 
used in Fig. H(a) but placed off axis for r = 0.2 and lattice 
size to = a = 0.01. The distribution of B particles is shown 
at instants t = 0, 0.2, 0.4 and 0.8. No B particle is in the wake 
of the cylinder after a time 1.0. (b) The same as Fig. ^(j) just 
in a more elongated frame. Note the downstream widening of 
the filaments. 

In fact, owing to a special symmetry, which is not 



present in the case of general obstacles, the flow is re- 
flection symmetric with respect to the x axis after a time 
shift of one-half. Therefore, the product distribution is 
of period 1/2 ||. 




FIG. 5. (a) The unstable manifold of the chaotic saddle in 
the reaction free flow generated by distributing 20.000 pas- 
sive point particles (black dots) on short segments along the 
local unstable direction of three basic fixed points of the La- 
grangean dynamics, and iterating them forward in time over 
several periods, (b) and (c) shows the B particle distribution 
at t = 10.0 just before and just after the reaction, respectively. 
Note the accumulation of B along the unstable manifold and 
note its sudden broadening from (b) to (c). The parameters 
are as in Fig 2. 

The outcome of the dynamics depends strongly on the 
initial position of the seed particle. If the initial B droplet 
is off axis [as in Fig. 0(a)], it does not penetrate the mix- 
ing region in the wake of the cylinder, and the initial 
droplet is just simply stretched before the whole amount 
of B is washed downstream. One can observe that the 
size of the compact patch B increases due to the autocat- 
alytic process as time goes on. Note that in this case no 
material B remains in the mixing region and the reaction 
dynamics dies out in any fixed observation region of finite 
size. This behavior is to be compared with cases where 
the droplet penetrates the mixing region. To sharpen 
the contrast, in Fig. |](b) we display the B distribution 
of Fig. |^(j) in a much longer region downstream. This 
clearly indicates that material B is now present at any 
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instant of time at any x-value in the wake. The gradual 
broadening of the stripes of product downstream is due 
to the autocatalytic feature of the process (and would 
not be present in the case of collisional reactions). 

In what follows we focus on such nontrivial cases in 
which the droplet penetrates the mixing region. To un- 
derstand the dynamics of Fig. ||, we recall that the tracer 
dynamics is governed by a chaotic saddle in the wake of 
the cylinder. Passive tracers coming close to the chaotic 
saddle spend a long time in the mixing region before be- 
ing advected away along the unstable manifold of the 
chaotic saddle [cf. Fig. g(a)]. Thus tracers having spent 
long time in the mixing region accumulate on the un- 
stable manifold. A comparison of Figs. |(i-j) and |(a) 
provides numerical evidence for the accumulation of ma- 
terial B in stripes of finite widths along this manifold. 

In order to gain more insight into the reaction dy- 
namics, Figs. ||(b) and (c) show the reagent distribution 
just before and just after the autocatalytic reaction takes 
place, respectively, in the steady state. In the first case, 
the B distribution has a rather scanty appearance, while 
right after the autocatalytic reaction most of the fila- 
ments of the manifold are washed out due to a sudden 
widening. The two pictures correspond to two differ- 
ent coverages of the fractal manifold. Just before the 
reaction, the unstable manifold is covered with stripes 
of average width £*(t), while just after the reaction with 
e* (0) > e* (t) . The sudden increase of the coverage width 
at certain times is due to our modeling of the chemical 
reaction as a "kicked" process. In the case of time con- 
tinuous reaction obtained in the limit r — > this feature 
is not present, but the fact that material B occupies a 
fattened-up fractal remains unaltered. We do not want 
to restrict our attention to the time continuous case, be- 
cause the more general setting might also be of relevance 
to biological phenomena, as discussed in the Introduc- 
tion. 



FIG. 6. Dependence of the number A* B {r)/a 2 of B parti- 
cles just before reaction in the steady state on the time lag 
t (all other parameters are as in Fig. 2, with the exception 
of the reaction range which is a = eo = 0.003, and of that 
a boundary layer of width 0.01 has been cut out around the 
cylinder surface). Note the sudden decrease around r x ~ 3 
(v x ~ 1-1), leading to a complete disappearance of B from 
the mixing region, indicating an emptying transition, around 

Tcrit — 6. 

One of the most interesting quantities to follow is the 
change of the number of B particles with the time lag r 
(or v) in the steady state, as shown in Fig. ^. Observe the 
monotonic decrease and observe that for relatively large 
r values (t > r cr i t ) no particle remains in the wake. This 
indicates the existence of an "emptying transition" . For 
reactions taking place rather seldomly, the effect of the 
advection by the background flow is so strong that no 
B particles survive the time lag r in the wake of the 
cylinder, and therefore, no reaction takes place in the 
mixing region. Such emptying transition occurs if r is on 
the order of the chaotic lifetime, and hence v cr n w 1. 



B. Collisional reaction: A + B — > 2C 

Initially the flow consists of material A, into which we 
inject reagent B continuously, along a line segment of 
length I perpendicular to the background flow in front 
of the cylinder. As time evolves, material C is produced. 
Figure U\ displays typical snapshots of the surroundings of 
the cylinder. The narrow stripes of constituent C (black) 
separating the areas occupied by background material A 
(white) and reagent B (grey) are clearly visible. Both B 
and C stripes are now pulled along the unstable mani- 
fold of the chaotic saddle forming lobes behind the cylin- 
der. This implies again that the reaction essentially takes 
place along this manifold [cf. Fig. ||(a)]. 
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FIG. 7. Time evolution for the collisional reaction 
A + B —* 2C. Snapshots (a-d) of the surroundings of the 
cylinder are taken at times 0.5, 1, 1.5 and 2, respectively, right 
before the reactions. The region —3 < x < 5, —1.5 < y < 1.5 
is shown. Initially the flow consists of material A, into which 
we introduce reagent B along the line segment x = —2, 
—0.05 < y < 0.05 (I = 0.1). White area denotes the back- 
ground material A, grey is reagent B, while black is the prod- 
uct C separating the former constituents. The reaction range 
is a — 0.011, £o = 0.01, while the time lag is r = 0.05. These 
imply an injection rate of 14.000 particles per unit time. 

Figure || shows the time evolution of the number of 
cells occupied by B and C. After a short transient, a 
saturation is reached. This means, that the number of C 
cells born in the reaction is the same as the number of C 
cells escaping due to the emptying dynamics. 
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FIG. 8. Time evolution of the number of cells occupied by 
B (shown by diamonds) and C (crosses) right before the reac- 
tion in the numerical experiment of Fig. 7. Note the asymp- 
totic stationarity reached by both components. 

As in the autocatalytic reaction, the approached sta- 
tionary value of C depends on the reaction range a and 
the reaction time v. The number of B cells also converges 
to a stationary value after sufficiently long time that de- 
pends not only on a and v, but also on the velocity v of 
the background flow, and on the width I since the amount 
of inflowing B per unit time is vqI. The number of B cells 
increases, in spite of the reaction that consumes B, due 
to the flow that brings initially more and more B parti- 
cles into the mixing region. A periodic pulsation can be 
observed after saturation has set in, just like in the auto- 
catalytic case. By increasing the time lag r, we observe 
that material B is not arranged along a fractal set in the 
wake, or does not reach the wake at all, due to the finite 
resolution of the grid. This already happens for not too 
large values of the grid size Eq . Because we are interested 



in the pronounced fractal structures formed, we keep the 
dimensionless reaction time considerable below unity in 
order to ensure that activity extends into the wake along 
the unstable manifold. 



IV. THEORY 

A. Autocatalytic reaction: A + B — > 2B 

Basic dynamics 

Let ,4b (i) = Ag(s) denote the area occupied by 
reagent B. Here s € [0,r] is the time after the nth reac- 
tion. Thus the total physical time is t = tit + s. During 
the time interval of length r only the dynamics of the 
chaotic flow controls the motion of the two components. 
In a fixed region of observation, say in a rectangle in the 
wake, Ag' \s) decreases according to the escape rate re of 
the chaotic saddle as 

^1^ = - K A%\ 0<s<t, (2) 
as 

provided the material is distributed along sufficiently nar- 
row stripes around the unstable manifold at s = 0. Solv- 
ing this equation we find that the area occupied by ma- 
terial B at the end of the interval r is 

A<g ) (T)=A<g ) (0)e-» (3) 

with v — kt [Q. 
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FIG. 9. Schematic diagram of the components' distribu- 
tion along the unstable manifold (bold line) before and after 
the nth reaction. The borderline between the components is 
denoted by solid (dashed) lines before (after) the reaction, 
(a) Auto-catalytic reaction, (b) Collisional reaction. Here all 
reagents A and B within a distance a are transformed into C 
in the course of reaction. The broadening of the half width of 
each C stripe at reaction n + 1 is denoted by > 0. 

In order to determine the new area ^ n+1) (0) right af- 
ter the (n + l)th reaction [which is expected to be larger 

than Ag(T)], we recall that after sufficiently long time 
from the onset of the reaction [when the exponential de- 
cay law (||) holds] the area of B is pulled into narrow 
stripes of more or less constant widths along the unsta- 
ble manifold [cf. Fig. @(a)]. In other words, the fractal 
unstable manifold is fattened up with material B which 
has a nonzero area. Let e(t) = e^ n \s) <C 1 denote the av- 
erage width of the stripes at time t = nr + s. This means 
that on scales longer than or equal to e^ n \s) the terri- 
tory occupied by B appears to be a fractal of the same 
dimension Dq as the unstable manifold. Covering the 
full area occupied by B at any time instant by squares 
of linear size e > e' n '(s), the number of boxes needed 
for this coverage behaves as N(e) = H.e~ D °. Here H 
is a constant characterizing the geometry, the so-called 
Hausdorff measure (or area) of the manifold. It can be 
obtained by determining the (unstable) manifold of the 
advection dynamics of the reaction-free flow. In what 
follows we assume H. to be known. 

It is worth noting that although the dimension is in- 
dependent on the instant at which the snapshot is taken, 
the Hausdorff measure is not. Since the flow is periodic, 
7i is periodic with the period of the flow: Tt(t) = H(t+1). 
(In fact, due the special reflection symmetry explained 
in Sec. III. A, H is periodic with period 1/2.) For conve- 
nience, we choose the period of the flow to be a multiple 
or a divisor of the time lag: Mr = 1, where M or 1/M 
is an integer, respectively. For time lags shorter than the 
flow's period, the period contains an integer number of 
reactions, otherwise the time lag is an integer multiple 



of the period. Thus H {n) (T) = H {n+1) (Q) = H {n) is M 
or 1/M periodic as a function of n. Since eW(s) is the 
smallest box size with which the fractality of the reagent 
can be felt, the area A^\s) of B can be written at any 
time as 



A { g\s) =H {n) (s) e< B )(«) 



2- Do 



(4) 



To determine the dynamics of the covering width 
eW(s) we observe that, at any reaction, there is a change 
due to the sudden increase of the product area [see 
Fig. 9(a)]. Since the filaments are locally smooth lines, 
this widening is proportional to the reaction range a. We 
can thus write 



(5) 



If the widening were exactly orthogonal to the manifold, 
and all the stripes occupied by B were nonoverlapping, 
the coefficient d^ would be 2. After some time, however, 
initially distinct stripes start to overlap. This leads to 
a change in the effective free surface available for the 
reaction. The phenomenological factor d' n ) introduced 
in Eq. (]^) describes both the effect of the geometrical 
shapes (stripes being tilted) and the effect of overlap, on 
the average. It depends on the time instant n, and on 
the flow parameters. The time evolution of the shape 
factor d^ n ' is unknown a priori. We shall see, however, 
that the theory becomes consistent with the numerical 
observations if, after a long time, the shape factor takes 
over the period of the flow. Thus we assume that d^ is 
also M or 1/M periodic in n. 

By taking into account Eqs. (|[^), the area A^ +1 \t) 
before the (n + l)st reaction can be written as 



A% +1) (T) = H in+1 ) 



,(n+l) 



(r) 



2- Do 



= e~ u n {n) [ £ (n+1) (0) 



2 -Dn 



From this we find that the average widths e( n+1 > (r) just 
before reaction (n + 1) is proportional to e^™ +1 '(0), the 
width right after the nth reaction: 



r («+i)( T ) = e (»+i)(o) 



l/(2-D ) 



Then from (||), a closed recursion relation follows for the 
widths just before the reactions: 



e (n+l)( T ) = 



l/(2-D ) 



(t) + o-d^] . (6) 



In view of Eq. (||) , this implies a recursion for the B area 
as 



A% +1) (T) = e-U A 



2 -Do 



(7) 
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This reaction equation is a discrete dynamics for A^ (r) 
expressing the amount of B before a reaction in terms 
of the amount before the previous reaction. Since the 
quantities and H {n) appear in a specific combination 
for this map, we have introduced the shorthand notation 



ee d {n) H {n) 1/{2 ~ Do) . 



(8) 



In contrast to the width dynamics (^J), the area dynamics 
contains all geometrical contributions via g^ , which is 
called therefore the geometrical factor. 

Equation (Q) belongs to the same class of dynamics 
as recursion relations of the type x n +\ — f{x n ), like 
e.g. the famous logistic map |3(|. It is one-dimensional 
and strongly dissipative, therefore it describes the con- 
vergence towards asymptotic motions which are repre- 
sented by attractors. If TiS n ^ = 7i and g^ = g is 
constant, a fixed point of the system is found from 
4*(r) = A { b +1) (t) = A ( £\t) in the form of 



A* b (t) =He* (2 - Da) = 



o-g 



,v/(2-D ) _ 1 



2- Do 



(9) 



This is the area occupied by reagent B right before a 
chemical reaction takes place in the stationary state. The 
area of B right after the reaction is a factor e lJ larger. 
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FIG. 10. Scaling properties in the autocatalytic reaction, 
(a) The area A* B = A* b (t) occupied by B in the steady state 
scales as A* B ~ a' 2 ~ D ° at different values of the time lag r 
(V t = 1, '*' r = 2, '+' t = 3) with Do = 1.64, in good 
agreement with the fractal dimension of the unstable manifold 
of the reaction free flow. Several runs were carried out with 
different values of the reaction range a, which is set equal 
to the lattice constant a — so- (b) Number N(e) of boxes 
of size e, covering the B material distribution right before 
the reaction, shows a fractal scaling. On the log 10 e 2 N(e) vs. 
log 10 e plot an exponent 0.4 appears implying Do = 1.6. The 
crossover value is e* ~ 0.003. Below this scale the fractal 
structure is completely washed out and a trivial exponent 
emerges. The parameters are a = 0.001, t — 1. (c) The small 
to medium time lag range of Fig. 6, (0 < r < 3) fitted to (^) 
with g = 5.5, k = 0.36, Do = 1.61, and a = 0.003. 

In the more general case when d^ n \ 7i^ n \ and thus 
g( n \ are periodic with some integer M, a limit cycle of 
period M is the attractor. The active process becomes 
thus synchronized to the underlying flow dynamics. Due 
to the linearity of Eq. (^|), and the fact that the factor 
e -vl(i-D a ) < ^ e \[ m [t C y C i e attractor of (Q) is always 



stable (in contrast to the logistic map). 

Equation (^) can be used to determine approximately 
the mean value of the geometrical factor jW and its oscil- 
lation amplitude about the mean in a temporally periodic 
steady state. For the case of Fig. 3, when the attractor 
is of period-5, we obtain g = 16 ± 3. Note that for in- 
teger values of r > 1 (when 1/M is integer, too), 
and WW are constant since the periodicity of the flow is 
unity. Then (S) holds again exactly. 

We have carried out a series of numerical experiments 
for such cases to carefully check the validity of (|]). 
Fig. |T(](a) shows the scaling of ln^4^(r) vs. lner. A scal- 
ing exponent of 0.36 emerges that is close to the theoret- 
ically predicted value of 2 — D ~ 0.39. 

A quantitative measure of the fractality of the B pat- 
tern can be obtained by computing the "box-counting" 
dimension. The number of boxes of linear size e that con- 
tain material B is expected to scale as e~ D °. Fig. |j"o|(b) 
shows a fractal scaling with 1.6 for e > e* and a sud- 
den crossover at e* due to the fattening-up of the frac- 
tal at small scales. This e* is the average width of the 
coverage of the unstable manifold and is found to be 
e* = 0.003 ± 0.0002. A comparison with (g) yields that 
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the shape factor is on the order of d = 0.6. 

Fig. |Io](c) shows a fit of Eq. @ to the measured values 
of A* B (r) . A good quality fit has been obtained with the 
geometrical factor g — 5.5 up to moderate values of the 
time lag r. 

Observe from formula ^ that, in the limit as r — > oo, 
the area of material B vanishes, A* b {t) — > 0. This does 
not necessarily mean that there is no material B left in 
the mixing region. In fact, there is always some material 
B left in this region, in the form of a fractal set of zero 
area. The r — > oo limit corresponds to the case where 
the B set literally coincides with the unstable manifold of 
the chaotic set. Since during a reaction the B set fattens 
up with average width proportional to a, the fattened-up 
B set should not have vanishing area as r — > oo. In fact, 
assuming that is constant, and using Eqs. (^) and 
(^), the limit r — > oo gives: 

lim A* B (0) = lim e"A* B {r) = (agf- D ° ^ 0. (10) 

r — >oo r — >oo 

Thus, taking the limit r — > oo in (^|) is only consistent 
with the chemical framework, i.e., with a continuous dis- 
tribution of particles. Of course, the numerics can deal 
with finite grid size only. However, when the grid size is 
small enough, we get a good agreement with our formu- 
las, for not too large time lags, as shown in Fig. [Tc|(c). 

Emptying transition 

When deriving the theory, we do not make use of the 
finiteness of the grid, which can also be interpreted as 
the finiteness of the particle size. It is this finite size ef- 
fect which leads to the emptying transition for large time 
lags. If during the advection dynamics of duration r the 
number of escaping particles exceeds the number of new 
B particles created at the next reaction, then the balance 
favors the total extinction of material B from the wake. 
This happens at a critical value for the dimensionlcss 
reaction-time, v cr it ■ 

We derive an expression for v C rit y i a defining the pro- 
ductivity of the reaction. The productivity of the chem- 
ical reaction in the steady state can be characterized by 
the ratio of newly-born to parent particles B as 



S=[A%(0)-A%(r)}/A B (T) = e"-l. 



(11) 



The productivity S, however, cannot be arbitrarily large. 
An absolute maximum S max exists, since the number of 
cells inside the reaction range a is limited. This implies 
that 



p v — 1 < <? 

For v smaller than 

Vcrit = In (1 + Smax) 



(12) 



(13) 



the productivity of the chemical reaction grows exponen- 
tially with increasing reaction time r or v. For v > v C rit, 



however, S does not grow further. Using Eq. ( O ) and 
the fact that S < S max , we have 

A* B (0)<(l + S max )A* B (T). (14) 

Inserting expression (||), this leads to 

A%(0)<(l + S max )A* B (0)e- v . (15) 

For v > v cr it the quantity (1 + S max )e~ v is less than 
1, therefore (15) can only hold for ^(0) = 0. Thus, 
for v > Vcrit, the area of reagent B quickly drops to 
zero in time, and remains only the background mate- 
rial A in the system in the steady state. In our case 
Smax = 8, thus the critical reaction time lag becomes 
Tcrit = 21n3/ft = 6.10. Indeed, this is confirmed by the 
measurements exhibited in Fig. |(| 

The definition of the productivity leads us to another 
important time value, the crossover time lag, v x ■ It rep- 
resents the threshold time lag of productivity for which 
the effects of finite particle size (and grid) come into play. 
Below this value, the filaments are 'tightly' covered by 
the B cells, so on average one particle is responsible for 
at most two newborns, S = 2. At v x , the coverage of the 



filaments is just about to break up, thus: 
is x = In (1 + 5) =ln3 , 



(16) 



which leads to r x = ln3/ft = T cr n/2 = 3.05. For values 
of t larger than r x , the theory described in the previous 
paragraph breaks down since it does not make use of the 
finite particle size. Indeed, a deviation from (^]) shows 
up in Fig. H by r = 3. 

The continuous time limit 

Before taking the continuous time limit t — ► 0, it is 
worth rewriting the reaction equation (^) in a different 
form. If the coverage of the manifold by B is relatively 
wide, i.e^, e'")(t) g^cr, we can expand the right hand 
side of (0) to first order to obtain 



A% +1 Ht)=A%Ht) e~» 



+ agW (2-D )e- v A%>(t) 



-0 



where 



Do-1 

2-Dn 



(17) 



(18) 



is a nonnegative number (it takes on zero only in non- 
chaotic flows) . This equation is equivalent to saying that 
the reaction gain is proportional to the perimeter C of the 
B area times the reaction range a |?j , and that this gain 
can be expressed to be proportional to a power — f3 of the 
area. [With the identification c = g i - n \2 - D ), Eq. @ 
coincides with Eq. (3) of Ref. 0.] Because exponent is 
negative, the equation contains a singular expression of 
the B area. The fixed point (0) for e* is, in general, of 
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the order of ga, and the condition leading to this form is 
typically not valid after a long time. Therefore, Eq. jl7| ) 
can only hold for a transient period before coming close 
to the attractor. If, however, the time lag is small enough 
and decreases, v <C 1, the region of validity of ([l7]) is in- 
creasingly longer. It is thus not surprising that in the 
continuous time limit, the differential equation obtained 
for the B area is the analog of the singular map (17): 



Ar 



-kAi 



g(t) (2 - D )v r (A B ) 



(19) 



Here A B (t) denotes -4^ /t) (0) or A% ,t \t) for r -> 0, and 

a 



v r = lim — 

T^O T 



(20) 



represents the velocity of the reaction front. 

A comparison with the case of a uniform flow [Eq. ([!])] 
shows that Eq. ( |l9| ) represents a novel form of reaction 
equation containing a negative power of the material con- 
tent due to the fractality of the unstable manifold. Its 
special case obtained for Dq = 1 describes a surface reac- 
tion in the presence of a single isolated hyperbolic orbit 
in the wake 



Ab = —kAi 



git) v r 



(21) 



Note that this is already structurally similar to the tra- 
ditional equation ([!]). But even for a single periodic orbit 
the decay is slower than without its existence because the 
time scale set by the escape rate (which is now just the 
Lyapunov exponent of the orbit) is typically much longer 
than the characteristic time of the flow. 

Altogether, Eq. (|l^) deviates in both terms from the 
traditional reaction dynamics (Q). The decay is slower 
while productivity is much faster than in an unstructured 
flow. Both effects are due to the presence of a chaotic 
saddle which produces an escape rate k smaller than vq 
and a dimension D bigger than one. 

Our findings concerning the new form of the reaction 
equation are similar in spirit to those of Muzzio and Ot- 
tino H . They considered the effect of filamentation on 
chemical reactions in closed containers without finding an 
explicit form for the reaction equations. The difference 
with us is due to the fact that we are studying now open 
flows with fractal patterns of Dq < 2, and our exponent 
j3 is therefore unique. 

The area occupied by B in a stationary state [g(t) = g] 
follows from Ab = as 



A* R = He* 



*2-D _ fff (2 - -Do) v r 



2-D 



(22) 



in accordance with the r — > limit of Eq. (||). Equa- 
tion (^||) has an important consequence for the velocity 
of a reaction with visible fractal properties. The latter 
can only be seen if e* is much less than unity (the cylinder 
radius) which implies, for g factors of the order of one, 
that v r <C K < vq. Fractal product distributions can only 



be expected for reactions which are slow on the lifetime of 
chaos. In this continuous time limit the emptying tran- 
sition does not occur, since e v — > 1, and therefore (|lj) is 
fulfilled. 

It is worth mentioning briefly that different continuous 
time limits are also possible depending on the assump- 
tions made on the model parameters. One option is to 
split the reaction range a into a 'deterministic' part a 
and a fluctuating 8a, i.e., a — a + 5cr. The reaction front 
velocity is then defined via <r, and 5a /t goes over as a 
noise term £ in the continuous time limit. So the reaction 
equation becomes a stochastic differential equation 



At 



-kA b + g{t) (2 - D ) (v r + (Ab) 



(23) 



Note that the noise term £ appears in this Langevin equa- 
tion in a multiplicative form. It's effect is thus also influ- 
enced by the fractality of the underlying manifold. 



B. Collisional reaction: A + B 



2C 



Basic dynamics 

The area occupied by the materials B and C at time 
s after the nth chemical reaction is denoted by A^ (s) 

and A^is), respectively. During the time interval r 
there is no reaction, thus only the chaotic advection and 
the injection of B governs the dynamics of the area of B 
and C as 



dA 



(«) 



= VqI — kA^\ 



dA, 



(n) 

c 



ds 



-kA 



in) 



(24) 



Here i>o denotes the velocity of the inflow far upstream 
from the cylinder, and I is the width of the stripe along 
which material B is being injected in the inflow region. 
The rest of the inflow consists of A only. After solving 
these equations we find the areas occupied by B and C 
right before the next reaction to be: 



A ( b\t) -^ l) (0)e- l/ + — (1 



vol 



(25) 



The reaction A + B — > 2C takes place if the distance 
between A and B is less then a. Thus, the amount 
A ( b +1) {0) of material B and _4? +1) (0) of C right after 
the (n + l)th reaction becomes smaller than Ag (r) and 

larger than A^ (r) , respectively. 

The geometry is now somewhat more involved than 
in the autocatalytic case. The branches of the unsta- 
ble manifold are covered with material B in stripes of 
average widths £^g\s). Adjacent to this, there are two 
stripes of equal widths e^ l \s) containing only material 
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C [Fig. 0(b)], while material A is outside. Thus, both 
materials B and C lie along a fattened-up copy of the 
unstable manifold. Note, however, that the amount of 
fattening- up is different for these materials. 

In an analogous way as for the autocatalytic reaction, 
we introduce again a phenomenological shape factor . 
Due to the effects of overlaps and non-straight geome- 
try of the unstable manifold, we assume that a reaction 

takes place at time (n + 1)t if Eq ' (t) is smaller than the 
range ad^ n \ Here dy 1 ' is again dependent on the time 
instant, being M or 1/M periodic in n after sufficiently 
long times. The broadening of the half width of each C 
stripe (which is the same as the corresponding decrease 
of the B stripes) in the (n + l)th reaction is 



.(«) 



(r). 



(26) 



This leads to the following change of the covering stripe 
widths: 



and 



4' l+1) (0)=e^ ) (T)-2C ( " +1) , 



4" +1) (0) =4? ) (t) + 2C ( " +1) 



(27) 



(28) 



Note that Eq. (|2^) is meaningful only if the difference on 
the right hand side is non-negative. This depends on the 
relation between the variables I and a. For a fixed I, a can 
be increased independently, and at a certain value, the 
increase in the width of the C-stripes exceeds the width 
of the B stripes during a reaction. From now on, we 
shall take a small enough, and assume that the situation 
in Fig. ^(b) is valid at all times during the process. 

Since the fractal scaling still holds for the actual width 
of B or C, at an arbitrary time, the areas are 



A { c\s) = 2H {n \s) 



» 



2-D 



_(n) 
'C 



2- Do 



(29) 



If the first equation is also applied to the inflow region, 
we have to assume that the width I of the injection of B is 
on the same order as £g(s), so that its two-dimensional 
character cannot yet be seen on this scale. Substitut- 
ing these into (p5|), we find two relations connecting the 
widths taken at different times: 



■4 B+1) (i-) = « (n+1) 



_("+!)/ 



e- u n in) 



.(n+1) 
~-B 



2- Do 



2-Do 



Vol 



(0) + 1^1(1 _ e -), 

K 



and 



A ( c +1 \t) = 2H {n+1) f4" +1) ( T ) 



2-D 



_(n+l) 



(0) 



2-D 



with H (n) = HW(t) = H^+^iO). In view of @ and 
(p8|) we obtain a coupled set of recursions for either the 

stripe widths or for (r) and A^ (t) . The dynamics 
of the covering stripe width ec for material C has the 
simpler form: 



(r) 



l/(2-_D ) 



2ad^-s^\r) 



(30) 



Consequently, the change of the half width is: 



(n+2) 



^(n+1) 



(n+1) 



c 

7^(n+l) 



l/(2-D ) 



(31) 



The recursion for the stripe width e # covering reagent B 
is then obtained from Eqs. (p5|), (p7|) and (|23) as 



4" +1) (r) 



4 n) (r)-2C 



(n+1) 



2-Dn 



l/(2-D„) 



(32) 



Finally, we find that the recursion relations for the 
areas are given by the following reaction equations: 



A in + l) {r) = Vjj L{l _ e _ v) + e _ v 



2 ^Oil) _ 2g (n) a 



2-D 



fA^(T)\ 1/i2 ' DoY 

2g (n) a _ *CM 



2- Do 



(33) 



Here denotes the same geometrical factor (||) as de- 
fined for the autocatalytic process. This form implies 
that the area dynamics for both B and C only depends 
parametrically on g^ n \ while the Eb and ec stripe dy- 
namics contain d^ and H' n ' as independent parameters. 

Note that the C-reaction (and ec) is decoupled from 
B (and £b), and component B simply follows the C- 
dynamics: the second of Eqs. ([33]) is independent of the 
B component, which may seem surprising at first sight. It 
was derived with the assumption that a is small enough, 
and that the average width of material B is large enough 
to furnish enough reagent for each reaction event: the B 
stripe is not totally consumed during an instantaneous 
reaction. Thus the dynamics of the reaction product C 
depends only on the actual width of the C stripes which 
separate A from B. Between two consecutive reactions, 
the average C width decreases, making possible a widen- 
ing at the next reaction, and so on. Thus, the presence 
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of B is necessary for producing C, but if the B-area is 
wide enough, the C-reaction becomes independent of B. 

In the case when the Hausdorff measure Ti and 
the shape factor d (and geometrical factor g) are n- 
independent, we find a stationary state. The average 
width of the C stripes can be given explicitly: 



s c (r)=e c (0) e-»l&- D °) 
Consequently, from Eq. ( |3l| ) 
C* = a d tanh 



2nd 



e v/{2-D ) + i ' 



2(2 -D ) 



(34) 



(35) 



follows. For the average width of the B stripes an implicit 
equation is obtained 



2-D Q 



vol 
kH 



(1 



[e* B {T)-2C?- Da . (36) 



Since in this equation all terms are expected to be of the 
same order of magnitude, we find that 



vol 
kH 



l/(2-D ) 



(37) 



holds. It means that the coverage width of the inflow of 
reagent B is of the order of the average coverage width 
in the mixing range, thus our earlier assumption on the 
validity of Eq. (| 2 9 [ ) is fulfilled. The fixed point of recur- 
sion ( B0|) and (|3l|)is an attractor for any parameter since 
e -u/{h-=b ) < t 

The fixed point expressions for the areas occupied by 
B and C are obtained as A* B {s) = He* B (s) 2 ~ Do and 
A* c {s) = 2He* c (s) 2 ' D ° , respectively. Thus, 



and 



2He 



* (2-D„) _ 

c — 



v l, 



2ag 



,u/(2-D ) _|_ i 



2- Do 



(38) 



^(rK + -(l~e") 

K 



l/(2-D ) 



A B (r) 



2ga 



A* c {r) 



l/(2-D„) 



l/(2-D ) 



(39) 



Note that although (|3^) is in some sense the analogue 
of (^), the particular ^-dependence is different: we have 
now a Fermi-Dirac distribution, while (^) was of Bose- 
Einstein type. This difference is due to the fact that the 
increase of ec in a reaction step is not a constant, rather 
it is proportional to £ which also depends on ec itself [cf. 
(psj)]. From these two relations a ^-independent form 
follows 



A% 

n 

1 _ ,W(2--Do) 



l/(2-D ) 



A b (t) 



^( T )\V(3-A.) 



l/(2-£>o) 



(40) 



which does not contain any free parameters. 

When the periodicity of H, d and g is pronounced, 
recursions (|||) typically possess a limit cycle attractor 
corresponding again to a full synchronization to the flow 
dynamics. 

The consistency of expressions (|38| ) , (|3^) and (^) with 
the numerics is verified in different ways (we fixed vq = 14 
and I = 0.1). The comparison of the results shown for a 
period 20 limit cycle steady state in Fig. 8 with the ex- 
pression ( |38| ) yields g — 42 ±4 for the geometrical factor. 
This value is also consistent with ([59|) and the measured 
B area. In simulations using different grid sizes £o, we 
find similar geometrical factor values that are slightly in- 
creasing with £o. We also evaluate the ratio of the left 
and right hand sides of the fit-free relation (40) at dif- 
ferent grid sizes £o = 0.007, 0.01 and 0.013 and time 
lags r = 1/80,1/40,1/20 and 1/10. The ratio is in all 
cases between 0.9 and 1.0 with an average 0.96 ± 0.05, 
0.95 ± 0.04 and 0.94 ± 0.03 for the e values investigated, 
respectively. In view of the fact that is not a con- 
stant, since the time lag is definitely below unity, the 
agreement, within an accuracy of 10 percent is satisfac- 
tory because this is exactly the amplitude of fluctuations 
in g. 

The continuous time limit 

Before taking the time continuous limit r — > 0, it is 
worth again rewriting the reaction equation ( |33| ) in a dif- 
ferent form. If the coverage of the manifold by B and C 
is large with respect to the amount of broadening, i.e., if 
M,£^(t) ;§> (( n+1 \ we can expand the right hand 
sides of ( p3| ) to first order in Q using ( |30| ) and ( |32| ) to 
obtain 
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(41) 



Here /3 is the same exponen t (18| ) as in the autocatalytic 
reaction. The fixed point (J34 ) for e* c is, in general, on 
the order of and the condition leading to this form is 
typically not valid after a long time. Therefore, Eq. (|4l| ) 
can only hold for a transient period before coming close to 
the attractor. If, however, the time lag is small enough 
and decreases, v <C 1, the region of validity of ( pl| ) is 
increasingly longer. 

Next we take the continuous time limit. Let fit) de- 
note the t — ► limit of any function which takes the 
form f( n \s) at time t = nr + s in the discrete time rep- 
resentation, with t kept fixed. Let us first consider the 
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recursion ( J3l] ) of ( and divide it by the shape factor d. 
The basic observation is that due to the minus sign in 
the first term on the right hand side, the recursion in the 
limit r — > does not go over into a differential equation 
for £/<£ Rather it shows that this ratio goes to zero lin- 
early with the time lag. Therefore it is natural to define 
a reaction front velocity as 



v r (t) = lim 



-►o rd(") 



(42) 



It should be emphasized that this property means that 
the broadening of the C stripes should tend to zero 
in the continuous time limit. Therefore, at any time 
£ c{t) = cd(t) holds. Thus, in contrast to the autocat- 
alytic reaction, the reaction range needs not to go to zero 
now, since another quantity, the actual broadening has 
to. In fact, er is finite, and is a measure of the front ve- 
locity. A substitution of this form into (|l|) yields the 
explicit expression 



(2 -D a 



dt^ 



(43) 



Since g is periodic with the flow's periodicity, the front 
velocity is also periodic. 

From recursion ( pd| ) we obtain for the differential equa- 
tions of the B and C areas in the limit r — > 



Ab = -kAb + v l - q(t) {AbY 
Ac/2 = - K A c /2 + q{t) {Ac/2) 



-a 



(44) 



where the coupling constant in both equations is the same 
q(t) = 2(2 — D ) v r (t)g{t). The two dynamics become 
entirely decoupled in this limit. Notice the singular terms 
on the right hand sides again. 

The fixed point value of the C area [for d(t) = d, 
H(t) = H, q(t) = q] is 



A* c = 2(q/n) 



2- Do 



(45) 



The B area has also a fixed point solution, but it is im- 
plicit: 



A*b = T~Ls\ 
The front velocity 



2-D 



2(2 -D ) 2(2-D )g 



(46) 



(47) 



is now a constant entirely determined by the reaction free 
flow's parameters and the reaction range. This implies 
ejr , = da in accordance with the v — > limit of equation 



V. DISCUSSION 

First we discuss a few aspects of the nonlinear dynam- 
ics which has not been mentioned in the text. 

The stable manifold is a fractal curve leading particles 
to the chaotic saddle. Points lying close to it reach the 
saddle rapidly, and their long lifetime is due to a resi- 
dence time in the neighborhood of the saddle. In other 
words, initial conditions close to or far away from the 
stable manifold lead to considerable lifetime differences 
only after the saddle region has been reached. Thus, al- 
though the role of stable and unstable manifolds seems 
to be symmetric in the problem, the residence about the 
stable manifold is not long enough for an enhancement 
in activity. Therefore, no fattening-up takes place along 
this manifold. Anyhow, a nontrivial reaction can only 
occur if the B (and C) initial conditions intersect the 
stable manifold. Otherwise, cases like the one shown in 
Fig. 4(a) occur without a fractal product distribution. 

The advection dynamics is known to have a consider- 
able nonhyperbolic component consisting of points lying 
very close to the cylinder surface. The nonhyperbolic 
part is characterized by a nonexponential decay (k = 0) 
and space-filling fractality (Do — 2) (^9|. Previous stud- 
ies [^,^lj have shown, however, that in the von Karman 
flow model the resolution allowed in computer simula- 
tions is still too crude to observe the nonhyperbolic ef- 
fects away from (but close to) the cylinder surface. The 
relative strength of the hyperbolic component ensures 
that on the time and length scales used we are able to 
work with a nontrivial fractal dimension and a finite es- 
cape rate. Apart from the nonhyperbolicity seen in the 
boundary layer around the obstacle, other nonhyperbolic 
structures are expected to be seen on very small scales 
only. The presence of a finite coverage with widths £*, 
however, prevents us to reach these scales. This is why 
the properties of the hyperbolic component play essential 
role in the full process. Therefore, in any fixed frame in 
the wake not overlapping with the boundary layer, the 
results of the theory presented here are expected to hold 

Next we summarize those features of our model which 
are believed to be general for active processes accompa- 
nied by weak diffusion in open flows with velocities faster 
than that of the reaction. 

• Active processes take place about the unstable mani- 
fold of the passive dynamics' invariant set. If the dynam- 
ics is chaotic, the manifold is a fractal and, consequently, 
the reaction leads to fractal patterns. 

• Although the fractal itself is a set of measure zero, 
the chemical products are of finite amount due to the 
fattening up process. 

• The fractal backbone results in an increase of the 
active surface, it acts as a catalyst, and generates an 
enhancement in activity as compared to flows with non- 
chaotic particle dynamics. 

• The inclusion of weak molecular diffusion in the 



14 



model can be regarded as a random walk superimposed 
on the deterministic advection. This would result in a 
further fattening- up of the fractal set, in a renormaliza- 
tion of a and thus of e*. 

• As a formal consequence of the fractal backbone, the 
product of the reaction obeys a singular scaling law. 

• In spite of the passive tracers' Hamiltonian dynamics, 
the active processes' equations are of dissipative charac- 
ter possessing attractors. 

• Most typically, a kind of steady state sets in after 
sufficiently long times, a state which is synchronized with 
the flow's temporal behavior. 

• Fractality is independent on whether the traditional 
reaction equations are linear or nonlinear since it is a con- 
sequence of the advection dynamics' strong nonlinearity. 

• The essential parameters for the chemical reaction 
in the flow depend on the parameters of the reactionless 
dynamics: the escape rate k and fractal dimension Do. 
These in turn depend on parameters (like the Reynolds 
number) of the underlying hydrodynamics. 

• The derivation of the reaction equations is simi- 
lar to the derivation of macroscopic transport equations 
from microscopic molecular dynamics. It seems that the 
presence of ever refining fractal structures (which can- 
not be observed directly with finite resolution) generate 
new terms in the reaction equation leading to observable 
macroscopic effects based on the fractal microstructures. 
They appear not only in the averages but also in moments 
if a stochastic description is used. 

All these features are expected to be present in realis- 
tic chemically or biologically active environmental flows 
observed on finite time scales. 
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